Irreducible tensor representation of MEG signals: Theory and applications 
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1. Introduction 

The MEG signal is a measurement of the magnetic 
flux through one or more coils, placed outside the 
head, in source free space. The magnetic flux is 
simply an integral of the magnetic field B over the 
directed area of the coil(s). Ideal MEG 
measurements (with no instrumental noise) must 
therefore obey the free space Maxwell’s equations, 
which in turn admit a special but well defined form 
for the magnetic field. More specifically, the 
magnetic field can be derived from a magnetic 
scalar potential U, which satisfies Laplace’s 
equations and can therefore be expressed as a series 
expansion of the well known spherical harmonic 
functions. Physics therefore defines the natural 
framework for the measurements, which in the 
quasi-stationary current source regime, allows 
separable treatment of spatial and temporal 
properties, with many convenient properties. We 
show how standard methods can be used to extract 
the (low order) expansion coefficients from any 
signal set and hence recast the measurements into 
the device independent U. We then show how this 
newly obtained U can be transformed into any real 
or virtual measurement set, or measurement related 
quantity. The resulting mathematical formulation is 
elegant and powerful, while the computations reduce 
finally to simple matrix manipulations. Some of the 
benefits of the methodology are discussed with a 
few of examples drawn from the many possible 
applications. 

2 Methods 

2.1 Basic concepts 

The MEG sensors are placed outside the source 
space (brain) and for the frequency domain of 
interest for neuromagnetic phenomena the quasi¬ 
static approximation holds. In this case, the 
magnetic field B satisfies V x B = 0, and hence can 
be expressed as the derivative of a scalar potential, 
B = VU, which satisfies Laplace’s equation 

v 2 u= o. 


2.2 Expansion series 

The space and time dependence of the scalar 
potential U can be written in separable form: 

U(f, t) = £ Yn ( r ) • c„ (0 = f(r) ■ c(t ), (2.2.1) 

n =0 

in terms of a set of time-independent functions: 

{y n : R 3 —>C; neN}. 

The magnetic field itself is by definition: 

B(f,t) = T(f)-c(t), (2.2.2) 

where 

r in {r) = -p 0 ^-y n {r), 

dr, 

r=r i -e,; i = 1,2,3; neN 
The matrix T is time-independent and contains 
information about the experimental design but it 
depends on the choice of expansion basis. Because 
of its time-independence, T needs be computed 
only once for a given sensor arrangement. The 
vector of coefficients c carries the complete 
temporal information of the experiment at hand, and 
is detector independent. The separable form of the 
expansion coefficients in the spatial and temporal 
domains can be exploited in many different ways. 

2.3 Determination of the expansion 
coefficients 

The measurement element associated with the 
sensor in is the magnetic flux, i.e. the integral of the 
magnetic field over the sensor area A m , which has a 
directed normal along d a m . In matrix notation the 
measurement set is given by: 

Hr,t) = Sl(r)-c(t), (2.3.1) 

with the elements of the time-independent and 
complex transfer matrix ft given by: 

ft m „ =-Po |d<i m ■ Vy„(r). (2.3.2) 

To get the vector of expansion coefficients c we 
need to solve the highly underdetermined, coupled 
linear system of equations 2.3. There are only M 



channels corresponding to the MEG device and an 
infinite number of coefficients. So we have to 
introduce an upper cutoff in the expansion N <M . 
The choice of N will have an influence on the 
accuracy of the expansion. If it is chosen too small, 
then details in the measurements will be unresolved. 
In general, allowing N = M is a bad choice because 
as the measurements increase, the additional 
information relates to very fine details, which sooner 
or later sink below the background instrumental 
noise. Setting N too high would therefore swamp 
the expansion with large components generated by 
amplification of relatively small noise components. 
There are a number of ways of choosing a cutoff 
with optimal properties, e.g. c = ft + -(j ), with 
ft + the Moore-Penrose pseudo-inverse of ft . Thus 
c can be determined directly from the measured 
signal (j) . Any feasible measurement or observable 
Zcan in turn be expressed as the product of the 
vector of expansion coefficients c and a static 
matrix M, containing both the device information 
and the observable characteristics: 

Z(t) = Mc(t)=MSl^(t). 


2.4 Spherical harmonics expansion 


The magnetic scalar potential is a solution of the 
Laplace’s equation and therefore expressible in 
term s of its eigenfunctions, the spherical harmonics. 
In the source-free space U satisfies Laplace's 
equation and hence has an irreducible representation 
in terms of the spherical harmonic functions, 
Y, m (6,(p); l = 0,...,L ; m = -l, as: 


n=l 2 +l + m 

where (p,9,(p) are the spherical coordinates of the 
measurement point. With this basic set of functions, 
the expansion in equation 2.2.1 can be rewritten as: 

U(pA<P) = ^^KP~ (M] -Yrm) , (2.4.2) 

and hence the elements of the transfer matrix ft in 
2.3.3 are given by: 

«„,(*») =~Mo Jd a,y\p~ u+1) ■ 7," (ftp)] . (2.4.3) 


2.5 Applicability to EEG 

A scalar potential V obeying Laplace’s equation is 
also directly related to the EEG signal, and therefore 
the framework provided by the spherical harmonic 
expansion is directly applicable to both EEG and 
MEG. For the case of EEG the measurements are 
taken on the surface of the scalp and in practice the 
shape and details of its conductivity distribution 
(which may not be easily determined) will influence 
the measurements. So, although in principle the 
mathematical formulation is applicable, in practice a 
completely model independent description may not 
be so easily available as for MEG. This is reflected 
by the earlier EEG studies using spherical harmonic 
expansion, where either a distribution on a spherical 
surface [1] or an explicit model for the conductivity 
[2] was assumed. We will restrict our examples to 
the MEG case only. 

3 Applications and results 

3.1 Virtual experiments and sensors 

A MEG measurement is frequently interrupted to 
allow the subject to relax. On other occasions, the 
same experiment is repeated on a different day. In 
these and many other situations the sensor geometry 
is (slightly) changed with respect to the head, 
making a direct comparison of the measured signal 
impossible. The same problem, but in a more acute 
form is encountered when one wishes to compare 
the results of the same experiment with the same 
subject but performed on different devices, perhaps 
with different numbers and types of channels. The 
formalism described earlier allows a neat solution to 
such problems, by expressing the signal that would 
have been measured with a detector B in terms of 
the signal that has been actually measured with a 
detector A : 

<j> B (t) = n B ft+ •£,(*). (2.5.1) 

This allows us to compare the results of a virtual 
experiment using the B -device (the real experiment 
used the A -device) with any real B -experiment. 
There is an analogous representation from classical 
optics: in the classical polarization experiment, an 
incident wave described by the field vector 
passes through a polarization filter characterized by 
ft^ index, the Polarizer, and afterwards through a 
second polarization filter characterized by ft B 
index, the Analyzer. The core concept is again the 



transformation of the device dependent signal (j> A 
(incident wave) via the fl) (Polarizer) into a device 
independent representation c. With that we can 
describe the signal (j) B (emergent wave) in any 
(virtual) experimental setup. 

Frequently, a weakly identifiable feature can be 
amplified by a linear combination of channels, 
making a virtual sensor. Such virtual sensors have 
been used with auditory and somato-sensory 
responses [3], Usually the peak at a particular 
latency in the average is used to define a linear 
signal combination, which amplifies the response 
from an area of interest. This virtual sensor can then 
be used to follow the responses in single trials. If 
such a response needs to be followed over long 
periods of time or on different days then one would 
like to have exactly the same linear combination of 
sensors, which demands that the sensors be placed in 
exactly the same position, which is rather difficult. It 
is however very easy to recover the same virtual 


sensor if we first transform the new measurements 
to the ones that would have been obtained if we 
placed the sensors at the old positions. 

3.2 Noise identification and elimination 

The spherical harmonics expansion offers a natural 
noise elimination mechanism for high spatial 
frequencies, because the spatial frequency content 
increases monotonically with the order of the 
expansion. Commonly, in MEG experiments the 
signal measured in an accurate channel is 
sufficiently smooth and the signal reconstruction is 
reasonably good, even for low expansion orders. In 
contrast, the signal measured in a noisy channel is 
very rough and the defect of the reconstructed signal 
will be significantly large at this sensor. A direct 
comparison between the actual and predicted signal 
provides an easy identification of noisy channels, as 
the examples in Figure 1 and Figure 2 demonstrates. 




Figure 1. The result for reconstruction of a low-noise channel (channel MRT11) on the left and a noisy 
channel (channel MRT12) on the right. The data were taken from an auditory experiment with the 151 sensor, 
whole-head Omega probe. The dotted (continuous) line shows the actual (reconstructed) signal. All 149 low- 
noise channels (including MRT11) were used for reconstruction. 
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Figure 2. The result for reconstruction for a noise-free channel, channel MRT11, and two noisy channels, 
channel MRT12, and MRT23, respectively. The MEG data were taken from the same experiment as for Figure 
1, but using the 148 low noise low-noise channels, which remain after the reconstructed channels are 
excluded. As before, the dotted and continuous lines represent the actual and reconstructed MEG signal, 
respectively. 


4 Discussion 

We have explored the general framework provided 
by the irreducible tensor representation of the 
magnetic field outside the head and used it to design 
a set of tools for manipulating the MEG signal. 
These tools are efficient because they allow the 
decorrelation of spatial features and dynamic 
aspects of the signal. The characteristic time- 
independent ft matrix needs only be computed 
once per experiment. The rest of the analysis 
reduces to a multiplication of a time-dependent 
vector with a static matrix. The method is also 
general because no assumption is made whatever 
about the generators or the conductivity model. The 
only assumption made is that all the sensors 
employed are outside the head, a fairly safe 
assumption to make. 

There are many theoretical as well as practical 
applications of the proposed tools. We have 
discussed just a few examples, which are of 
considerable practical importance, namely the 


identification of noisy channels, or noisy segments 
in the MEG signal and the construction of virtual 
sensors and virtual experiments. 
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